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Motivated by recent experimental reports, we carry out a Fermi liquid many-body calculation 
of the interaction induced renormalization of the spin susceptibility and effective mass in realis- 
tic two dimensional (2D) electron systems as a function of carrier density using the leading-order 
'ladder-bubble' expansion in the dynamically screened Coulomb interaction. Using realistic material 
parameters for various semiconductor-based 2D systems, we find reasonable quantitative agreement 
with recent experimental susceptibility and effective mass measurements. We point out a number 
of open questions regarding quantitative aspects of the comparison between theory and experiment 
in low-density 2D electron systems. 
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I. INTRODUCTION 

It is well-known that the mutual Coulomb interaction 
between electrons could cause substantial quantitative 
modification of thermodynamic properties (e.g. effective 
mass, specific heat, compressibility, magnetic susceptibil- 
ity) in an interacting electron liquid. This is the so-called 
many-body renormalization of the Fermi liquid parame- 
ters, which has been studied extensively in three dimen- 
sional metals^ and in two dimensional (2D) semiconduc- 
tor structures^ for a very long time. At zero temperature 
(or more generally at low temperatures, T/Tp <^ 1 where 
Tp = Ep/ks is the Fermi temperature) the many-body 
Fermi liquid renormalization for a quantum electronic 
system is entirely determined by the electron density (n) 
with the dimensionless density parameter rg being de- 
fined as the average inter-electron separation measured 
in the units of Bohr radius: Vs = (7rn)~^/^/aB where 
qb = nh^ /{me^) is the effective Bohr radius for a back- 
ground (lattice) dielectric constant k and a band mass m 
- this definition of applies to 2D (in 3D: r^ cx n~^/^) 
with n being the 2D electron density. It is easy to see 
that our Tg-parameter is proportional to the ratio of the 
average Coulomb potential energy (i.e. the interaction 
energy) to the noninteracting kinetic energy, and as such 
the system is strongly interacting at large (low den- 
sity) and weakly interacting at small (high density). 
We emphasize that our definition of does not depend 
on the spin (or valley degeneracy in 2D) of the system. 
Studying (and comparing theory with experiment) the 
density dependence of various Fermi liquid parameters 
in interacting electron liquids has been one of the most 
important and active many-body research areas in con- 
densed matter physics in 3D metallic systems (as well 
as in normal IIe-3, a quintessential Fermi liquid albeit 
with a short-range inter-particle interaction) and more 
recently, in 2D electron systems confined in semiconduc- 
tor structures. The 2D systems have the distinct advan- 
tage of the density being a tunable parameter so that the 
density dependence of the Fermi liquid renormalization 



can be studied directly. In this paper, we theoretically 
consider the density-dependent many-body renormaliza- 
tion of the 2D electronic spin susceptibility and effective 
mass, a subject of considerable recent experimental activ- 
it}i2i^i^iSiL2i2iiSiii in a number of different semiconductor 
heterostructures with confined 2D electron systems. 

There has been a number of experimental papers 
appearing in the recent literature reporting the low- 
temperature (< lOOmK) measurement of the susceptibil- 
j^j,3j4j5j6j7j8j9 ^^j^^j effective massiS^ii in 2D electron sys- 
tems as a function of carrier density in the Ts w 1 — 10 
parameter range. Although some aspects of the data in 
different experiments (and more importantly, the inter- 
pretation of the data) have been controversial - most 
especially on the issue of whether there is a spontaneous 
density-driven ferromagnetic spin polarization transition 
at low carrier densities in 2D systems - the experimen- 
tal reports convincingly establish a strong enhancement 
in both the spin susceptibility and effective mass as a 
function of decreasing (increasing) carrier density n (in- 
teraction parameter Ts). This strong enhancement of the 
susceptibility with decreasing carrier density has been 
demonstrated in 2D electron systems confined in (100) Si 
inversion layers, in GaAs heterostructures, and in AlAs 
quantum wells. The typical enhancement in the low tem- 
perature susceptibility is in the range of a factor of 1 — 4 
for Ts ~ 1 — 10. In addition to the strong low-density 
enhancement of the measured low-temperature suscepti- 
bility, there are several other interesting and intriguing 
features in the experiments. The susceptibility enhance- 
ment shows a modest dependence on the induced spin 
polarization (or equivalently, a change in the spin de- 
generacy) in the 2D system - since the susceptibility is 
typically measured^^ by applying an external magnetic 
field to produce a spin splitting in the 2D system, the 
susceptibility is invariably measured in the presence of 
finite spin polarization (and the extrapolation to zero 
spin polarization may not always be uniquely reliable). 
Another interesting observed recent feature, reported in 
AlAs quantum well systems, is that the susceptibility en- 
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hancement is larger in the single-valley semiconductor 
system rather than in the multi-valley system in con- 
trast to the expectation based on just exchange energy 
considerations (since the multi-valley system is naively 
expected to be more "dilute" as the electrons are be- 
ing shared among different valleys). The enhancement 
of effective mass with decreasing carrier density was ob- 
served in Si inversion layer and GaAs heterostructures. 
This measurement is performed by examining the tem- 
perature dependence of the low-temperature Shubnikov 
de Haas oscillation magnitudes. The effective mass mea- 
surement turns out to be more difficult in general. Very 
recently, a detailed and extremely careful experimental 
measurement of density dependent effective mass, using 
several internal consistency checks, has appeared in the 
literature^i. 

In this paper, we provide an excellent qualitative 
and reasonable quantitative (realistic) theoretical under- 
standing of the density-dependent 2D spin susceptibil- 
ity and effective mass measurements at low temperature. 
Our work, in contrast to much of the existing theoreti- 
cal literature on the topic, fully incorporates the realistic 
quasi-2D nature of the electron systems (i.e. the fact that 
these systems have finite widths in the transverse direc- 
tion normal to the 2D plane of confinement, which consid- 
erably modify the Coulomb interaction between the elec- 
trons) which is of substantial quantitative importance in 
the experimental parameter regime. We also study the 
spin- and the valley-degeneracy dependence of the cal- 
culated 2D susceptibility, obtaining in the process qual- 
itative agreement with the recent experimental finding 
on the anomalous valley dependence of the susceptibil- 
ity in 2D AlAs quantum well structures. (The spin, gs 
and the valley degeneracy, g^ dependence of the suscep- 
tibility enters through the 2D density of state, which is 
proportional to gsgv) Our best quantitative agreement 
with the experimentally measured 2D effective mass is 
obtained for the so-called "on-shell" self-energy approxi- 
mation. 

In Sec. ^we lay out the theory and formalism of our 
calculation. We present the calculated results for spin 
susceptibility in Sec. IIIII and effective mass in Sec. IIVI 
for different realistic 2D systems. In the end in Sec. El 
we discuss upon various issues related to our results and 
their relations with the experiments. 



II. FORMALISM 

We calculate the T = (paramagnetic) spin suscepti- 
bility X* and effective mass m* of a quasi-2D electron sys- 
tem by using the many-body perturbation theory tech- 
nique. It has been knownH for a long time that for 
an electron system interacting via the Coulomb interac- 
tion, the most important terms are associated with the 
long-range divergence (the so-called 'ring' or 'bubble' di- 
agrams^*^) of the Coulomb interaction, and as such an 
expansion in the dynamically screened Coulomb inter- 



action (with the screening implemented by the infinite 
series of 'bubble' polarization diagrams) is the appropri- 
ate theoretical framework. Such an expansion is, in fact, 
asymptotically exact in the — > high-density limit, 
and is known to work well qualitatively for > 1 al- 
though its precise regime of validity can not be deter- 
mined theoretically and has to be ascertained by com- 
paring with experiments. The theory becomes progres- 
sively quantitatively worse as increases, but unless a 
quantum phase transition intervenes, there is no specific 
Ts value which necessarily limits the qualitative validity 
of the theory. Motivated by the encouraging fact that 
this leading-order expansion in the dynamically screened 
interaction (involving an infinite resummation of all the 
bubble diagrams in the electron self-energy calculation) 
seems to provide a good quantitative description^^ for 
the thermodynamic properties of the 3D simple metals 
(rs « 3 — 6), we apply here the same theory for interact- 
ing electrons in quasi-2D semiconductor structures. 

(c) (d) 




(ii) (iii) 



FIG. 1: (Color online.) (a) The RPA ladder- bubble series for 
the interacting susceptibility with the bold straight line the 
interacting Green's function and the bold wavy line the dy- 
namically screened interaction; (b) the noninteracting suscep- 
tibility; (c) the Dyson's equation for the interacting Green's 
function in terms of the noninteracting Green's function and 
the self-energy; (d) the self-energy in the leading-order ex- 
pansion in the dynamical screening; (e) the Dyson's equation 
for the dynamically screened interaction in terms of the bare 
Coulomb interaction (thin wavy lines) and the polarization 
bubble; (f) a charge fluctuation diagram which does not con- 
tribute to spin susceptibility; (g) Landau's interaction func- 
tion. 

Our many-body diagrams for the interacting suscepti- 
bility is the so-called "ladder-bubble" series as shown in 
Fig-D^ this is a consistent conserving approximation for 
the susceptibility. Direct calculation of these diagrams 
turns out to be difficult for the long ranged Coulomb 
interaction. However, at T = 0, Landau showed that 
X* can be equivalently expressed through the Landau's 
interaction function /(k, k') (shown in Fig. ^g)) asi^: 

-.=—.+77^1 feme, (I) 

X m* \2ttY J 
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where x is the PauH spin susceptibiHty, fe{d) = /e(k, k') 
with k and k' on-sheU (i.e. /2m = k'^ /2m{= Ep)) 
is the exchange part of Landau's interaction function 
(Fig- Eg) (i)); S is the angle between k and k'. Similarly, 
the Landau theory expression for the effective mass m* 
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Note that the spin independent exchange Landau's inter- 
action function /e(k, k') is responsible for the difference 
between the ratio x/x* ^-nd m/m*~. 

An equivalent, and easier way to derive the effec- 
tive mass is through calculating quasiparticle self-energy 
and obtaining its momentum derivative. The self-energy 
within RPA can be written asi^ 
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dv 



_fq_ 
(27r)2 J 2711 e(q,j/) 



Go(q+k,i.+a;), (3) 



where Vq = F{q)2'Ke'^ / q is the bare electron-electron 
interaction and F{q) the quasi-2D form factor for the 
electron-electron interaction^ which will be described in 
detail later. Note that this form factor also appears in 
the expression of dynamical dielectric function e(q, ly) 
through Vq. As mentioned earlier, a significant quantita- 
tive feature of our theory is the inclusion of the realistic 
quasi-2D Coulomb interactions, in our calculation which 
substantially reduces (compared with the pure 2D case) 
the quantitative many-body renormalization effects. In 
Eq. 
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is the bare Green's function with np the Fermi distri- 
bution function and = k'^ / (2m) — Ep with Ep the 
non-interacting Fermi energy. It is showni-4 that the in- 
tegration along real axis in the expression of self-energy 
(Eq. (O) can be deformed onto imaginary axis, which 
avoids the singularities along the real axis and makes the 
integration easier. The contour deformation also breaks 
the expression of self-energy into separate terms that 
correspond respectively to contributions from the spin- 
dependent and spin-independent part of the Landau's 
interaction function shown in Fig. ^ and is very useful 
for us to derive the expression for susceptibility as shown 
below. The expression of the real part of the self-energy 
can then be written as 



Re E(k, ^) = - J (^^'90(2"^^ + 4 - |q - kp) 

(2^)2^'''%(q,^q_k-c.) 
e{2muj + kl-\q- kp) - e{kl - |q - kp) 



d^q 
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where kp is the Fermi momentum. The effective 
mass is derived from the expression of the real part 
of the quasiparticle self-energy by m/m* = 1 + 
{m/kp)-^Re I](k, ^k)|/c=feF^- Combining this with 
Eq. 1^, we have 

V 1 /" Tn d 

^ = l+(^ J UO)dO+--ReEiK^.)\,.,.. (6) 

It is not difficult to show that the second term 
of Eq. jSJ accounts for the contribution from the 
spin-independent exchange Landau's interaction func- 
tion /e(k,k') (Fig. Eg)(i)), and therefore the term 
(27r)^2/ fe{9)d9 in Eq. © exactly cancels the momen- 
tum derivative of the second term in the self-energy 
Eq. ([SJl. Hence the expression of x/x* only contains con- 
tributions from the k derivatives of the first and third 
term in Eq. (|SJ|. After converting all the expressions in 
terms of the dimensionless parameter r^, and using 2kp, 
AEp, 2m as the momentum, energy, and mass units, the 
expression for x*/x = 9*m*/{gm), where x*{x), 9*{9), 
m*{m) are respectively the interacting (noninteracting) 
spin susceptibility, the Landau g-factor, and the effective 
mass, is given in our theory as 
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where a = ^Jg^gg/A with g^ [gs) the valley (spin de- 
generacy); A — x"^ — x"^ — u,B = 2xu, R = \J + i?^; 
X — ql{2kp) and u — uj/{AEp)] e{x,iu) = 1 + 
arsgvgsF{x)[l I {2x) — A + R/{2^/'^x'^)\ is the imaginary 
frequency dielectric function. For 2D quantum well, the 
form factor F{q) is given by 
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where d is the width of the infinite square-well potential 
of the quasi-2D system. For heterostructure quasi-2D 
systems (e.g. Si MOSFETs), the form factor is 



^, , , Kins, 8 + 996+3^262 

^^''^-^^^—J 8(1 + ,6)3 
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'2(l + g6)6' 
(9) 

1 /3 

with b = [Kscfi^/(487rrn^e2n*)] defining the width of 
the quasi 2D electron gas, Ksc and Kins are the dielectric 
constants for the space charge layer and the insulator 
layer, is the band mass in the direction perpendicular 
to the quasi 2D layer, and n* = ridepi + with ridcpi 
the depletion layer charge density and n the 2D electron 
density. We choose ndepi to be zero in our calculations, 
since it is unknown in general (finite small value of ridcpi 
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do not change our results). From Eqs. |^ and ^ we 
have 

m X , otrs F{x) 
m* X Jo xe{x,()) 

To avoid any possible confusion, we emphasize that the 
definition of our spin susceptibility is the derivative of 
the magnetization with respect to the applied magnetic 
field, and therefore x* is well-defined even in a non-zero 
parallel magnetic field B. ft is important to note and 
easy to show that, within our approximation in which the 
spin-orbital effect can be neglected, the parallel field B 
dependence of the spin susceptibility x* ^-nd the (7-factor 
g* manifests only through the spin-degeneracy factor g^, 
which is present in the expression (|7|). 

Motivated by recent experimental studies, we have di- 
rectly evaluated the interacting susceptibility as a func- 
tion of density at different spin-((7s = 1,2) and valley- 
degeneracy {qt} = 1, 2 for AlAs quantum wells) for 
three different 2D semiconductor systems: n-Si(lOO) in- 
version layers; n-GaAs gated undoped heterostructures; 
modulation-doped AlAs quantum wells. We also calcu- 
lated effective mass for n-GaAs heterostructures moti- 
vated by a very recent experimental reportii. In the 
rest of this paper, we present and discuss our calculated 
results for x* Ix ^-^d m* /m in light of the recent experi- 
mental data in 2D semiconductor structures. 



III. RESULTS FOR SPIN SUSCEPTIBILITY 
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FIG. 2: (Color Online.) Calculated spin susceptibility of Si 
MOSFET ideal and quasi-2T> systems with Qs = 1 and 2, to 
be compared with experimental data from Refs. 00- 

In Fig. 13 we show our (100) Si MOSFET results (tak- 
ing a valley degeneracy — 2) for the calculated suscep- 
tibility comparing with the recent experimental results. 
We show two sets of results corresponding to the realistic 
quasi-2D system and the strict 2D system. In each case 



we show the calculated susceptibility for both the = 1 
(due to the presence of an external magnetic field) and 
the normal gs = 2 situation. Since the experiments are 
invariably carried out in the presence of finite external 
magnetic fields, the experimental results probably corre- 
spond to the region in between the gs — I and gs = 2 
theoretical curves (for the quasi-2D system). A more de- 
tailed discussion is made in Sec. on the reason why 
we are comparing the experimental data with our gs = 1 
and gs = 2 theoretical results. There are three important 
points we make about Fig. [21 (1) The quasi-2D results 
are lower than the 2D results by a factor of 1.5 to 3, 
and the relative difference is much larger at low carrier 
densities since the effective quasi-2D layer width is larger 
at lower 2D densities; (2) the theory gives a reasonable 
description of the experimental data, - in particular, the 
experimental data points lie very close to the region be- 
tween the theoretical gs = 1 and gs = 2 quasi-2D sus- 
ceptibility results; (3) the pure-2D results disagree with 
experiments. 
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FIG. 3: (Color online.) Calculated spin susceptibility of GaAs 
heterostructure ideal and quasi-2D systems with — 1 and 
2, to be compared with experimental data from Refs.Q 

In Fig. 121 we show our theoretical susceptibility results 
for electrons confined in quasi-2D GaAs heterostructure 
(valley degeneracy gv — 1) comparing with the recent 
measurements of Zhu et. ali. Again, the agreement 
between our results and the experimental measurements 
is very good for the reahstic quasi-2D calculations. Note 
that at the low carrier densities of interest in the GaAs 2D 
system, the quasi-2D effects are extremely strong (almost 
a factor of 5!), and the quantum Monte Carlo (QMC) cal- 
culationi^ with which the experimental results were com- 
pared in Ref . are completely inapplicable since they are 
for a strict 2D system rather than a quasi-2D system. In 
fact, any agreement between the measured susceptibility 
and the QMC results for the ideal zero-width 2D sys- 
tem should be taken as a spectacular quantitative failure 
for the QMC theory in the low-density regime. The rea- 
son for such a strong quasi-2D effect in GaAs is that at 
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very low densities (n ~ lO^cm"^), the transverse quasi- 
2D width of the electron wavefunction is extremely large 
(> 500A) so that the Coulomb interaction is substan- 
tially suppressed compared with the strict 2D limit re- 
sult. (This demonstrates that it is misleading to compare 
quasi-2D experimental results with the strict 2D theory 
as is often done in the literature!). 



r 




01 23456789 10 

n(10^^ cm-^) 



FIG. 4: (Color online.) Calculated spin susceptibility of AlAs 
quantum well ideal and quasi-2T> systems with (js = 1 and 2 
and Qv = 1 and 2, to be compared with experimental data 
from Refs.0. Effective k = 10. We use two different values of 
band mass: m = 0.21me and m — 0.46me. 

In Fig.0]we show our theoretical results for AlAs quan- 
tum wells. The AlAs quantum wells used in Ref. |30 are 
rather narrow (width ~ 50 A), and the additional self- 
consistent confining potential arising from the 2D elec- 
trons themselves produces even stronger confinement of 
the carriers thus further narrowing the effective quasi- 
2D width. Therefore, we show only the strict-2D results 
for the 2D AlAs system in Fig. 0] for the sake of clarity. 
Our corresponding quasi-2D results are somewhat below 
the theoretical curves for the strict-2D case. In Fig. ^ 
we show our theoretical results for two different valley 
degeneracy (gy = 1 and 2) situations. Our theoretical 
investigation of the valley degeneracy dependence of the 
2D spin susceptibility behavior is necessitated by the puz- 
zling recent experimental finding^ of an interesting val- 
ley degeneracy dependence in AlAs systems, namely, the 
many-body enhancement for x*/x is larger for the val- 
ley occupancy of 1 than the valley occupancy of 2. (Our 
results for Si MOSFETs and GaAs heterostructures in 
Figs. 1 and 2 are for valley degeneracy values g^, = 2 and 
1 respectively, consistent with the known band structure 
for Si and GaAs.) 

In Fig. ^ the results for m — 0.21me, which corre- 
sponds to the transverse AlAs band mass rrit, are in 
reasonable quantitative agreement with the experimen- 
tal AlAs quantum well results, including the anomalous 
finding of x* Ix having a stronger many-body Fermi liq- 
uid enhancement for the lower valley occupancy of 1 than 



for the higher occupancy of 2. We note that at high 
enough electron densities (> lO^^cm"^), this anomalous 
valley dependence of the spin susceptibility would disap- 
pear according to our theoretical calculations with x* Ix 
for gy = 2 being larger than that for = 1, consistent 
with one's naive expectations. It may be more appropri- 
ate to usei& a larger 2D band mass in our calculation 
given by m = 0.46me = y'mtmi with mi = rUe the 
longitudinal band mass. Use of this band mass results 
in larger theoretical values for spin susceptibility than 
that was measured in experiment, but the trend and the 
characteristic of valley-degeneracy dependence remains 
the same. The exact cause of the quantitative difference 
between our m = 0.46me results and the experiments 
probably lies in the yet unknown details of the samples 
and experimental procedures, and are therefore not clear 
to us right now. 



Why does the valley occupancy dependence of the 
many-body susceptibility enhancement act in such a 
way? The reason lies in many-body correlation effects be- 
yond the naive exchange energy considerations. Within 
an exchange-only Hartree-Fock theory, in fact, the 2- 
valley occupancy state would have a higher suscepti- 
bility enhancement than the 1-valley state. But corre- 
lation effects are important in the system, and at low 
enough carrier densities the 1-valley state turns out to 
have stronger many-body effects than the 2-valley state. 
Our theory, which is essentially a self-consistent dynam- 
ically screened Hartree-Fock theory, includes correlation 
effects demonstrating that at low enough carrier densi- 
ties, X* /x could be enhanced in the 1-valley state over 
the 2-valley state as has been experimentally observed. 
Another way to understand this is through the screen- 
ing effect. As gy increases, Fermi momentum decreases 
as fci^^ oc l/y^, which favors many-body renormaliza- 
tion. However, the screening effect increases with in- 
creasing gy since qtf oc which tends to decrease 
the renormalization effect. These two effects are com- 
peting with each other. At low densities, the screening 
effect is predominant, and therefore a smaller gy results 
in much less screening than a bigger gy and hence pro- 
duces larger susceptibility renormalization. On the other 
hand, at high densities the screening effect is less im- 
portant, and a smaller gy results in a larger Fermi en- 
ergy and a smaller renormalization effect. In particular, 
the valley degeneracy dependence of many-body effects 
should qualitatively correlate with the dimensionless pa- 
rameter qTF/(2kF) oc gv^"^: for large (small) qtf / {2kp)^ 
low (high) density, smaller (larger) values of gy produce 
larger renormalization. Note that the experimental re- 
sults in Ref. 9 do not show any strong gs dependence in 
contrast to our theoretical results in Fig. 01 We discuss 
this puzzle in Sec. Ivl below of our paper. 
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FIG. 5: (Color online.) Calculated effective mass of GaAs 
heterostructure ideal and quasi-2T> systems with Qs — 2 and 
(?„ = 1, to be compared with experimental data from Refs.lTll 



IV. RESULTS FOR EFFECTIVE MASS 

Although the main task of our current work is the cal- 
culation of interacting 2D susceptibility to compare with 
the recent experimental results2i^i^iSiL2i2ii£, a very recent 
experimental reporUi of 2D n-GaAs effective mass mea- 
surement allows us to apply our recently developed the- 
ory^'* for the quasiparticle effective mass enhancement to 
2D GaAs in order to compare with these experimental 
results. For details on the quasiparticle effective mass 
theory we refer to our recent publication^"*. 

In Fig.Elwe show our theoretical effective mass results 
for electrons confined in quasi-2T) GaAs heterostructure 
(valley degeneracy — \) comparing with the recent 
measurements of Tan et. al}^ . The agreement between 
our results and the experimental measurements is very 
good for the realistic quasi-2D calculation. Again the 
ideal-2D effective mass results turn out to be much larger 
than the experimental data, emphasizing once more the 
importance of the quasi-2D effect on the many-body 
renormalization of physical quantities in such systems. 
We emphasize that, while at high densities the quasi-2D 
results are rather close (with a 10 — 20% difference) to 
the pure 2D results, at low densities this difference could 
be as large as a factor 2. 



DISCUSSION 



effects through the standard variational approximation^ 
for heterostructures. We believe that our quasi-2D elec- 
tronic structure model is extremely reasonable, but some 
uncertainty and error probably arise (particularly at low 
carrier densities) from our lack of knowledge about the 
depletion charge density in these systems. 

A rather important factor in the comparison with the 
experimental data that is left out of our consideration 
is the effect of orbital couplingii of the in-plane com- 
ponent of the external magnetic field invariably present 
in the experimental measurement of the spin susceptibil- 
ity. At low carrier densities, when the quasi-2D width 
of the 2D layer is large, such a magneto-orbital coupling 
could have substantial effects on the 2D effective massif. 
This effect is, however, entirely of one-electron origin, 
and we assume, somewhat uncritically at this stage, that 
the magneto-orbital effect drops out of the susceptibility 
enhancement factor x* /x since the enhancement involves 
a ratio of the interacting and the non-interacting effective 
mass both of which will be affected in a similar manner 
by the magneto-orbital effect. In any case, the magneto- 
orbital effect is negligibly small for Si MOSFETs and 
AlAs quantum wells because of their tight quasi-2D con- 
finement, and in the n-GaAs structure the ratio x*/x 
should not be much affected by the magneto-orbital cou- 
pling. Only at rather low densities in GaAs heterostruc- 
tures the magneto-orbital effect may play a role. 
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FIG. 6; (Color online.) (A) Comparison between the results 
of spin susceptibility calculated using QMC method^- and our 
results calculated within RPA for a 2D electron system with 
Qa = 2,g-u = 1. Our results are for both ideal and quasi- 
2D systems as shown; the QMC results are for ideal 2D sys- 
tems. (B). Comparison between our RPA many-body theory 
and the QMC results— for quasiparticle effective mass. Re- 
sults are shown for strict 2D (both on-shell and off-shell RPA 
self-energy approximations as well as QMC) and quasi-2D 
(on-shell RPA self-energy approximation) cases. Our effec- 
tive mass results shown in Fig. 5 correspond to the on-shell 
approximation which is the consistent approximation for the 
Fermi liquid effective mass within the RPA self-energy ap- 
proximation. 



Before concluding we make some critical observations 
and comments on our theory and its implications and, 
more particularly, on the comparison with the recent ex- 
perimental results. One of our important conclusions is 
that the inclusion of quasi-2D form factor effects is essen- 
tial in understanding the 2D susceptibility results. As 
such, an important issue is the realistic nature of our 
model where we have included the quasi-2D form-factor 



Although our quasi-2D model is reasonably realistic 
and accurate, our many-body spin susceptibility calcu- 
lation is necessarily based on approximations since the 
problem of an interacting quantum Coulomb system can- 
not be solved exactly. We emphasize in this respect that 
QMC calculations, while being accurate as a matter of 
principle, turn out as a matter of practice to be often 
unreliable due to various approximations (e.g. node fix- 
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ing, back-flow correlations, etc.) and inaccuracies (e.g. 
finite size effects, numerical errors) inherent in the QMC 
technique. Also the QMC calculations leave out quasi- 
2D effects, and are therefore quantitatively incorrect. In 
Fig. I^lwe make a comparison between our result and the 
QMC resultif for susceptibility in an ideal-2D electron 
system with — 2 and — 1. In the intermediate 
and large region, the difference between the results of 
RPA and QMC methods turn to be very large. But we 
note that our quasi-2D results are quantitatively quite 
close to the QMC pure-2D results. (A similar situation 
seems to hold for the effective mass calculations also.) 
This demonstrates the care that is needed in comparing 
theory with experiments in this problem. 

In comparing theory and experiment for the suscepti- 
bility and the effective mass in 2D electron systems it is 
important to emphasize the fact that the experimental 
measurements for x* a-nd m* in 2D systems are typi- 
cally not thermodynamic measurements, but are essen- 
tially transport measurements of thermodynamic quan- 
tities in an applied magnetic field which spin-polarizes 
the system through the Zeeman splitting. For example, 
the 2D susceptibility is measurecl2iiiSiSiL2i2ii£ by either 
monitoring'' the full spin polarization of a 2D system in 
an applied parallel field or through the clever coincidence 
techniquei^ in a titled magnetic field. We want to raise 
three important issues related to these experimental tech- 
niques and their relevance to our theory. First, strictly 
speaking these techniques does not provide a true mea- 
surements of the spin susceptibility x* that is defined 
through x*{B) — dM/dB with M the magnetization 
and B the applied magnetic field. In all these techniques, 
the experimentally measured quantity is actually M/B, 
which we now call x*{B)- As has been explained in the 
experimental literature^i^i^iSiL^iS, these studies measure 
indirectly the occupancies of spin up and down levels 
through transport or Landau level coincidence measure- 
ments. The spin susceptibility is then inferred from these 
indirect spin polarization studies by using the simple non- 
interacting spin-dependent density formula: In the trans- 
port measurement technique, g* ^sBp — Ep where Bp 
is the so-called saturation (parallel) field for 'complete 
spin polarization'. In the Landau level coincidence tech- 
nique, g* hbB = itUc = iehB±/m* with i half integer 
or integer. In either case, the derived spin susceptibility 
is exactly x*{B) = M/B that we have just discussed. 
{B = Bp in the transport measurement technique, and 
B corresponds to the Landau level coincidence magnetic 
field in the Landau level coincidence technique.) It is 
important to note that x*{B) ^ X*{B) ^ X*{B = 0) 
for B >> 0, and these three quantities coincide in the 
B limit. The precise definition of spin susceptibility 
was never explicitly made clear in the previous experi- 
mental literature. Second, even though the experimen- 
tally observed x*{B) is not exactly the thermodynamic 
quantity x* = X*{B = 0), in normal circumstances they 
are close to each other. Especially, if the magnetization 
curve of the system is smooth (which normally it will 



be), x*{B) is in between x*{B = 0) and x*{B = Bp) 
which are the two extreme case we are considering in 
this work. As we mentioned before, since the magnetic 
field dependence of x* only manifests itself through the 
spin degeneracy factor, we can say that the experimen- 
tally measured x*{B) should be in between x*{9s — 2) 
and x*(ffs = l)j which is the justification for our com- 
parison of the experimental data with our gs = 1,2 the- 
oretical results. In Ref. 0, the issue of finite field and 
hence the spin degeneracy factor dependence of the sus- 
ceptibility is made more clear. However, in this case not 
only the partial spin polarization need to be considered, 
the effect of the finite perpendicular magnetic field and 
Landau levels on the susceptibility need to be further in- 
vestigated, which is beyond the scope of this work. Third, 
the difference between x*{B) and x* rnay be helpful in 
understanding some of the discrepancies between the ex- 
perimental data and our theoretical findings. For ex- 
ample the recent measurements^ in 2D AlAs quantum 
well systems find a strongly valley-degeneracy-dependent 
spin susceptibility which, however, demonstrates essen- 
tially no spin degeneracy dependence. Since the valley 
index is essentially a pseudo-spin index, theory would 
predict the same valley and spin degeneracy dependence 
of susceptibility unless there is strong spin-orbit or valley- 
orbit band structure effects in the system. Therefore, in 
the Ts region where the spin susceptibility has a large 
dependence on the valley degeneracy, the spin degener- 
acy should play an equally big role, as shown in Fig. ^ 
This theory-experimental discrepancy may very well be 
a result of the difference between x*{B) and x* ■ This 
factor might have also caused the quantitative discrep- 
ancy between our RPA results (using m = 0.46me) and 
the experimental results shown in Fig. ^ Other possible 
explanation could also lie in one-electron band structure 
physics, not in many-body theory. There have been no 
QMC calculations revealing the spin and valley degener- 
acy dependency of the spin susceptibility. 

The many-body approximation we use in our work, 
namely the renormalized interaction calculated within 
the infinite series of bubble diagrams and the susceptibil- 
ity calculated in the infinite series of ladder diagrams us- 
ing the appropriate dynamically screened interaction (i.e. 
the ladder-bubble approximation for x*)i corresponds to 
the leading-order expansion in the dynamically screened 
interaction. Such an expansion is asymptotically exact 
in the high-density limit, and is known to work well in 
the low-density limit as well. Unfortunately, its regime 
of quantitative validity is unknown, and it is likely to 
become progressively quantitatively inaccurate as in- 
creases. It should, however, be qualitatively valid for 
larger Ts values as long as the interacting electron system 
remains a Fermi liquid. Our reasonable agreement (with- 
out adjusting any parameters) with existing experiments 
shows that the theory remains well-valid in 2D-systems 
at least up-to ~ 7, which is consistent with the success 
of the corresponding 3D theory in the metallic densities 
{ts ^ 3 — 6). Since no systematic and uncontrolled many- 
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body theory approximation is available for susceptibility 
beyond the ladder-bubble series approximation carried 
out in this paper, we can only test the validity of our 
theory through the direct comparison with experiments, 
and on this ground the theory seems to be well-justified. 

Our quasiparticle effective mass results (Fig . are 
based on the theory developed by us in Ref. [ij recently. 
This theory is the RPA theory, based on the leading- 
order self-energy calculation in the dynamically screened 
Coulomb interaction. While the theory is asymptotically 
exact in the high density r — > limit, it is not based on 
an Ts expansion and is essentially a self-consistent field 
theory. Unfortunately, however, like many other self- 
consistent filed theories, the level of validity of this theory 
for large Ts (i.e. the strongly interacting regime) is un- 
known except that it is expected to remain qualitatively 
validi*i^*i^ for > 1 unless there is a quantum phase 
transition. Since we have discussed the approximation 
scheme and the validity of our RPA effective mass theory 
in some details in Ref. Q, we do not repeat those argu- 
ments here. We point out that the effective mass results 
given in Fig. [51 use the "on-shell" self-energy approxima- 
|-jQj,i3,i4^ which has been argued— to be the correct dy- 
namical approximation consistent with the Landau Fermi 
liquid theory as long as the self-energy is obtained in 
the leading-order dynamically screened interaction (i.e. 
RPA). Our "off-shell" approximation for the 2D effective 
mass (shown in Fig.|H{B)) differs from the corresponding 
"on-shell" results by a factor of 2 or more at low densi- 
ties, and show different trends as well. The fact that the 
"on-shell" effective mass results agree much better with 
experiment than the "off-shell" results is additional evi- 
dence in support of the "on-shell" approximation being 
the correct one within RPA self-energy scheme. One can 
try to "improve" upon RPA by including local field cor- 
rectionsiS to the dynamical electron polarizability (i.e. 
bare bubble of RPA) which, in some crude manner, sim- 
ulates the incorporation of higher-order vertex correc- 
tions in the theory. But such local field corrections are 
uncontrolled, and probably inconsistent, since many di- 
agrams in the same order are typically left out. We are 
therefore unconvinced that the inclusion of local field cor- 
rections in the theory is necessarily an improvement on 
RPA. The great conceptual advantage of RPA is that it 
is a well-defined approximation that is both highly phys- 
ically motivated (i.e. dynamical screening) and theoreti- 
cally exact in the high-density (r^ 0) limit. Attempted 
improvement upon RPA through the arbitrary inclusion 
of local field correction may neither be theoretically jus- 
tifiable nor more reliable. Keeping these caveats in mind, 
we mention that several previous worksMiiS have shown 
that local field corrections to the many-body properties 
of two dimensional electron gas turn to be very small in 
the density range of our current interest. 

In has recently been discussed in detail by us^"* that 
RPA is not necessarily a high-density theory although 
it is exact in the high density limit. The ring diagram 



approximation, which is at the heart of RPA, is an expan- 
sion in the dynamically screened interaction. Under some 
circumstance RPA can be very loosely thought o^i*i^ as 
an expansion in an effective parameter Ts/ {rs + C) where 
C > 1 is a constant. This implies that RPA could in 
some situation turn out to be decent approximation in 
the Ts > 1 case. Indeed 3D metals, with 3 < < 6, 
seem to be reasonably well-describedi by RPA. Leading- 
order vertex corrections to RPAiSiiiiiS typically produce 
only small quantitative corrections in the Ts > 1 regime, 
again empirically justifying the validity of RPA. There 
is obviously a lot of cancellation among the higher-order 
diagrams (in particular, between the self-energy and the 
vertex diagrams) leading to the good empirical success 
of RPA consisting only of the bare bubble diagrams. 
But this empirical success is physically well-motivated, 
arising from the long-range nature of the inter-electron 
Coulomb interaction with the decisive point of physics 
being the dynamical screening properties of a Coulomb 
liquid. In an interacting system with short-range bare in- 
teractions, e.g. neutral He-3, the ring-diagrams play no 
special role, and RPA is not a meaningful approximation. 
In a strongly interacting system, such as a 2D electron 
system with rg > 1, one must use a careful and critical 
comparison between theory and experiment as the prin- 
cipal guide in deciding the validity of any theory. By 
this standard RPA appears to be a reasonable approxi- 
mation even for > 1. There are many other theoretical 
techniques in condensed matter physics which work well 
beyond their putative regime of validity ~ for example, 
the dynamical mean field theory (DMFT), which is ex- 
act in infinite dimensions, seems to provide good results 
in 2 and 3 dimensional strongly correlated systems, and 
the local density approximation (LDA) , which is exact for 
very slowly varying density inhomogeneity, works well for 
the band structure of real solids with their rapidly vary- 
ing densities. Similar to DMFT and LDA, RPA works 
well because it is fundamentally a non-pcrturbative the- 
ory which accounts for some key ingredient of the inter- 
action physics, namely, screening in the situation under 
consideration in this paper. 

We conclude by emphasizing that the most significant 
implication of the excellent qualitative and quantitative 
agreement between our theory and 2D spin susceptibil- 
ity and effective mass measurements is that interacting 
2D electron systems, which have been of much interest 
recently, remain Fermi liquids down to reasonably low 
carrier densities {rg ^ 7 or so at least), and many-body 
perturbative techniques, coupled with a realistic quasi- 
2D description for the electronic structure, provide a very 
good qualitative and quantitative model for these sys- 
tems. There is no need to invoke any non-Fermi liquid 
concepts to explain the recent experimental results on 
quasiparticle renormalization effects^i^i^iSiLSiSiiSiii. 

This work is supported by the NSF, the DARPA, the 
US-ONR, and the LPS. 
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